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Abstract — The peak energy resolution of germanium detectors 
deteriorates with increasing neutron fluence. This is due to hole 
capture at neutron-created defects in the crystal which prevents 
the full energy of the gamma-ray from being recorded by the 
detector. A phenomenological model of coaxial HPGe detectors is 
developed that relies on a single, dimensionless parameter that is 
related to the probability for immediate trapping of a mobile hole 
in the damaged crystal. As this trap parameter is independent of 
detector dimensions and type, the model is useful for predicting 
energy resolution as a function of neutron fluence. 
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I. Introduction 



A gamma-ray traversing a germanium crystal loses energy 
mainly by the production of Compton electrons (or a pho- 
toelectron in the case of low-energy gamma-rays), which in 
turn lose energy by the production of electron-hole pairs at 
roughly 2.96 eV each. The electrons and holes are collected 
at the outer and inner electrodes covering the annular surfaces 
of the cylindrical crystal, thereby recording the energy of the 
gamma-ray. The full energy is not recovered, however, when 
holes are trapped at negatively charged defects in the crystal 
created by a flux of fast neutrons. Thus the energy resolution 
of gamma-ray peaks deteriorates with fast neutron fluence. 

In practice, the gamma-ray energy is determined by mea- 
suring the current flow, induced by the moving electron and 
hole charges, between the electrodes. In this phenomenological 
model, the gamma-ray energy is instead calculated from 
the integrated charge at the two electrodes, induced by the 
electron-hole pairs. 

This paper is laid out as follows. In section II the relation 
between an electron-hole pair and the induced charges at the 
electrodes is derived: this is essentially the concept of the de- 
tector In section III expressions for the spatial density of hole 
traps (regions of fast neutron damage), and for the interaction 
cross-section of those traps, are derived. In anticipation of the 
computer implementation of the model, section IV describes 

Manuscript received 2011. This work was supported by the Defense Threat 
Reduction Agency and the U. S. Department of Energy. It was performed at 
the Idaho National Laboratory, a DOE laboratory operated by Battelle Energy 
Alliance under DOE Idaho Operations Office Contract DE-AC07-05ID14517. 

The authors are with the Nuclear NonproUferation Division, Idaho National 
Laboratory, R O. Box 1625, Idaho Falls, ID 83415 USA. E-mail addresses: 
clinton.vansiclen@inl.gov, edward.seabury@inl.gov, jayson.wharton@inl.gov, 
gus.caffrey@inl.gov. 



the stochastic methods for choosing the location of electron- 
hole pair creation and subsequent location of the hole trapping. 
In section V the computer program, that combines these pieces 
of the model to produce spectral line peaks, is described. 
Finally, model results are compared to measurements, made 
in our laboratory, of the energy resolution of the 1332 keV 
^°Co spectral line obtained by a p-type detector irradiated by 
a fast neutron fluence up to 10^ cm^'^. 

The primary purpose of this work is to produce a predictive 
model of peak energy resolution as a function of fast neutron 
fluence. It is a phenomenological model, meaning that critical 
parameters are taken from experiment rather than calculated 
from first principles. The (dimensionless) critical parameters 
in this case are a^, which is related to the hole trap cross- 
section, and at, which is essentially the average number 
of hole traps created by a fast neutron as it collides with 
atoms in the germanium crystal. (These are both introduced 
in section III.) Values for these parameters are obtained by 
reproducing, with the model, the data of R. H. Pehl et al. 
lU. In that work, two HPGe coaxial detectors (one n-type, the 
other p-type), fabricated from the same crystal, were irradiated 
simultaneously with fast neutrons from an unmoderated ^^^Cf 
source. 

For the convenience of the reader, SI units (and derived 
units) are used in any calculations: electric potential (f) in volts 
(V); electric field E in V m^^; force in newtons (n); charge in 
coulombs (C). Note that 1 V = 1 n m C~^. Further, the product 
e (w V) = w eV. The electric charge e = 1.602 x 10^^^ C; the 



vacuum permittivity eo 



.854 X 10 
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Where the symbols ± and ^ occur, the top/bottom sign is 
used for p-type/n-type detectors. The coaxial detector is an 
annular cylinder so it is natural to use cylindrical coordinates 
(r, z): the radius of the inner electrode is i?o; the radius of the 
outer electrode is Ri ; the axial coordinate z = at the top of 
the coaxial detector (the end pointing at the radiation source), 
and z > Q within the detector 

II. Induced charge at the electrodes due to 

ELECTRON-HOLE PAIRS 

As an electron and hole have opposite charge, it is not until 
they move apart under the influence of the driving force gE(r) 
that charge is induced at the electrodes. The variable q is the 
charge of the particle (so equals +e/— e for holes/electrons) 
and E(r) is the electric field at the particle location r between 
the electrodes. In the case of p-type detectors, the outer contact 
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is positively biased so that the field E(r) is negative; thus 
holes move towards the inner contact at Rq and electrons move 
towards the outer contact at Ri. In the case of n-type detectors, 
the inner contact is positively biased so that the field E(r) 
is positive; thus holes/electrons move towards the outer/inner 
contact. 

What charges are induced at the electrodes by an electron- 
hole pair after separation? Expressions for these are easily 
found by use of Green's reciprocation theorem [2J: For a 
given arrangement of electrodes, if </> is the potential due to a 
volume-charge distribution p and a surface-charge distribution 
a, while (p' is the potential due to other charge distributions 
p' and a', then 



a'(l)da. (1) 



In the "primed" system, set p' — <t' = 0, so (j)'{r) is the 
solution to Laplace's equation in cylindrical coordinates with 
boundary conditions 0' = 1 at surface S, which is inner 
contact Rq, and </)' — at outer contact Ri. Then the right 
hand side equals zero, so the "unprimed" system (which is the 
annular cylinder with a point charge at r and a surface charge 
at the inner contact) obeys the relation 







/ a(f)' da ~ / p (f)d^x + 
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>(j)'d^x — ~ crda. 
'V J s 

Recognizing that p (the charged particle at r) is q times the 
delta function produces the relation q(p'{r) = — Jg crda. Thus 
the induced charge at the inner contact due to a charged 
particle at r is 



Qo{r) - ~q(t>'{r) 



Mr/Ro] \ ^ ln[fli/r] 
\n[Ri/RQ\) ^ln[i?i/i?o]' 

Similarly (but setting </)' = 1 at surface S, which is outer 
contact Ri, and 0' = at inner contact i?o), the induced 
charge Qi{r) at the outer contact due to a charged particle at 
r is 

ln[r/i?o] 



(4) 



'In [Ri/RoY 

Note that the induced charges Qo{r) and Qi{r) at the two 
electrodes are opposite in sign to the charge q of the particle 
at r, and that the sum Qo{r) + Qi{r) always equals —q. 

Now consider an electron-hole pair created at ri. As the two 
particles are of opposite sign and both reside at r^, no charge 
is induced at the two contacts. However, as the two mobile 
charges separate and move radially to and rh, respectively, 
the charge Qq is induced at the Rq contact and the charge Qi 
is induced at the Ri contact, where 

ln[r^/re] 



Qo = Qo{re) + Qo{rh) = e 



Qi = Qi{re) + Qiirh) 



In [Ri/Ro] 
In [rh/ve] 



(5) 



(6) 



'lTi[Ri/RoY 

It is these induced charges Qq Qi that account for 
the current between the two contacts due to an electron- 
hole pair (note that an incremental change SQo = —6Qi 



always, as electric charge flows from one contact to the other). 
Note that initially — = so the induced charges 
Qo — Qi =0, and that when both mobile charges successfully 
reach their respective contacts, the induced charges Qo ~ Te 
and Qi — ±e, as expected (top/bottom sign indicates p- 
type/n-type detector). As the current between electrodes is 
just the transfer of charge from one to the other, the result 
Qo| — \Qi\ = e allows the full creation energy of the 
electron-hole pair to be recorded. 

For a p-type detector, when the hole is trapped at but 
the electron successfully reaches Ri, the charge induced at the 
central contact is Qo = -^plllLHilLiA. 

at the outer contact is Qi — l. , ,„ ,„ , 

\n[Ri/Ro] 

case Qo = Qo {rh) and Qi = e + Qi (rh) < e, meaning 
that IQol = IQil < e so not all of the pair creation energy is 
recorded. For an n-type detector, when the hole is trapped at 
but the electron successfully reaches Rq, the induced charges 



. ,„ ,p 1 and the charge induced 

ln[Ri/Ra] = 

MMuA Note that in this 



are < 



, _ „ Hrh/Ro] 
- ^\n[Ri/Ro] 



and Qi — — e 



and again not all 



. _ . 'ln[fliARo] 

of the pair creation energy is recorded. By examining the two 
expressions for Qo (or Qi), it is evident that in general the 
average value (|Qo|) for a damaged p-type detector will be 
less than that for a damaged n-type detector, since a spatially 
uniform flux of gamma-rays will produce more electron-hole 
pairs near the outer contact than near the inner contact, so 
resulting in more hole trapping near the outer contact (thus 
{\n[Ri/rh]) < (In [r/j/i?o]))- This effect translates into less 
energy being attributed to an incident gamma-ray by a p-type 
detector than by an n-type detector. 

The induced charges (currents) at the contacts are related 
to the energy of the incident gamma-ray in a straightforward 
way. When the electron and hole reach their terminal locations 
at Te and r^, respectively, the induced charges Qo and Qi 
correspond to an energy recorded by the detector, where 
Q = IQol = IQil and e is the average energy needed to create 
an electron-hole pair (this is the energy needed to elevate an 
electron in the valence band into the conduction band). For 
germanium, e = 2.96 eV at 77 K An incident gamma- 
ray of energy that produces n electron-hole pairs will thus 
record an energy n-^e, where the average value (Q) is taken 
over all the n pairs. 

In the implementation of this model in a com- 
puter code, a value for n is chosen for each inci- 
dent gamma-ray from the Gaussian distribution p{n) — 



(2^2 



-1/2 



exp 



(n — riy) I (2it2) . The average value 



= E^jt [so for example, the 1332 keV Co gamma-ray 
produces (on average) = 4.5 x 10^ electron-hole pairs in 
a germanium detector]. The variance = F ■ where F is 
the Fano factor (and is approximately 0.13 for Ge detectors 
|3|). To obtain values for n from this distribution, it is most 
convenient to use the Box-MuUer method |4|: Consider the 
"standard normal" distribution — (27r) exp [—1/^/2] . 
Then the two random variables v\ and V2 will both have the 
standard normal distribution and will be independent, where 



V2 



(—2 Xmxx^^"^ cos (27rx2) 
(—2 Inxi)^^^ sin (2-kx2) 



(7) 
(8) 
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and xi and X2 are random numbers taken from the uniform 
distribution on (0, 1]. Then the desired value n ~ + av, 
where v is either v\ or vi calculated from Eq. Q or Eq. 
Note that the FWHM of the Gaussian distribution p(n) 

1/2 

is 2 (2(7^ In 2) . The FWHM of the corresponding energy 



peak is then 2 (2Fn^ In 2)''-' • e = 2 (2i^i;^eln2)''^ So for 
example, in the absence of neutron damage, the FWHM of 
the peak corresponding to the 1332 keV ^"Co spectral line is 
1.686 keV, while that of the peak corresponding to the 122.1 
keV ^^Co spectral line is 0.510 keV. 



III. Hole trapping 

The fate of mobile holes (and electrons) in the neutron- 
damaged crystal is determined by the spatial distribution of 
traps, and the trap cross-section. These two functions are 
derived in turn. 

The distribution of particle traps should be uniform over 
a z-slice (thickness dz), since the source of gamma-rays (at 
which the detector is pointed) also acts as the primary source 
of neutrons. What then is the trap density Pq(z) (the number 
of traps in the infinitesimal volume A - dz at axial position z, 
where A is the cross-sectional area of the crystal)? Consider 
that the mean free path of an incident neutron is /. That is, the 
probability that the neutron first collides a distance [z, z + 
from its entry point into the crystal is exp [— z//] • [dz/l). Then 
N exp [—z/Z] • [dz /I) is the number of neutrons that first collide 
a distance [z, z + dz] from their entry point at z = 0, where 
N is the number of neutrons incident on the crystal. Assuming 
that the collision produces a trap, the trap density 



Pg{z) = at ■ Nexp [~z/l] ■ {dz/l) / {A ■ dz) 

^ r in 
= at ■ -^exp [-z/l\ 



(9) 



where N /A is the neutron fluence, and the dimensionless 
parameter at is the (average) number of traps created by a 
fast neutron as it collides with atoms in the germanium crystal. 
Note that at should be somewhat larger than 1, to account — in 
a crude way — for any subsequent collisions by the neutron. 

What is a reasonable value for atl For a neutron incident 
on a detector of length (thickness) L, the probability that no 
collisions occur over the distance L is exp [—L/l], meaning 
that for a neutron fluence N /A, the fraction 1 — exp [—L/l] 
undergo at least one collision. L. S. Darken et al. after 
irradiating a 3 cm thick Ge crystal with neutrons, conclude "A 
fast neutron flux of 10^" cm^^ produces about 2 x 10^ cm"'^ 
disordered regions of various sizes". Setting / = 6 cm fE\ and 
L = 3 cm, the fraction of incident neutrons that underwent 
collisions in the crystal was 0.3935. Each of those neutrons 
was then responsible for (2 x 10^ x 3) / (0.3935 x lO^'') = 
1.5 traps per (collided) neutrons. Thus in general at ^ 1.5. 

As charged particles, holes and electrons may be trapped at 
defects with opposite charge. It is believed that hole traps are 
large disordered regions with large negative charge Q, and 
so have a large effective cross-section, while electron traps 
are much smaller (perhaps point defects) and so have a much 
smaller cross-section. In any event, the trap cross-section (jq 
must be roughly the size of the local distortion, due to the 



electric charge of the defect, of the applied electric field E(r). 
An expression for cr^ can be derived as follows. 

For simplicity, use 2D Cartesian coordinates, and place the 
defect (with charge Q) at the origin, and set the no-defect 
electric field E — Ex.. Then the potential (j)*{r) at the point 
r = (x, y) for this system is 4> {x) + ^ {x^ + V^) where 
{x) is the potential at r in the absence of the defect and e is 
the permittivity of germanium. This produces the electric field 



E* (r) = -V0* {r) = {E 



^12 
— X \X 

e 



y 



-3/2 



^ / 2 

—y [x 



y 



-3/2 , 



(10) 



The electric field lines in the absence of the defect are 
directed parallel to the x-axis; in the presence of the charged 
defect at the origin, they are bent towards the origin. Those 
field lines that terminate at the defect are particle paths that 
lead to trapping. Clearly the field lines are seriously bent 
towards the defect when the magnitude of the y-component 



Q 



-3/2 



of the field E* (r) exceeds that of the x- 



y[x-+ y 

component E + ^x (x^ + y^) that is, when x 



y" < 



Thus the trap cross-section <t ^ ny^ oc 



eE 



and 
As 



the electric field E(7-) in the crystal has a radial dependence, 
+ dz] o'g(r) = aq ■ ^E{r)~^ where aq is a dimensionless parameter. 

What are reasonable values for aql L. S. Darken et al. |8| 
estimate cross-sections ah ^ 10^^^ cm^ and Ue ^ 10~^^ cm^. 
A typical value for E{r) (the magnitude of the electric field 
in the crystal produced by the bias potential at an electrode) 
is 125 kV m-^ Thus a,, - 0.1 and a^ - 10~^. 

Due to the higher production of electron-hole pairs near 
the outer contact it is preferable to have a larger electric 
field E there as well to reduce the trap cross-sections aq. This 
shaping of the electric field is accomplished by doping p-type 
and n-type germanium detectors with (electron acceptor) boron 
and (electron donor) lithium, respectively. These dopants pro- 
duce an intrinsic space (free) charge density p = —ep^ in the 
case of p-type detectors and p = epo in the case of n-type 
detectors ||9l, where pa/d is the density of acceptor/donors. A 



typical value is pa/d — 10 cm^''. 

The potential (f){r) between the contacts satisfies Poisson's 
equation, V^(/) = — p/e, where e — 16eo is the permittivity of 
Ge. This equation 

p{r) 



(r-<j>(r) 
r dr \ dr 



(11) 



is solved for (j>{r) given the boundary conditions, which are the 
applied potentials at the outer and inner contacts. For p-type/n- 
type detectors, the outer/inner contact is positively biased. The 
electric field between the contacts is then E(r) = — V(/)(r). 
In the usual case that the charge density p has no radial 



dependence, (j){r) = cj) {Rq) 



4e 



Rl) +Aln[r/Ro] and 



E(r) = (^r — ^) r where the constant 

^^_ ^iR^)-q>iRo) + £{Rl-Rl) 
ln[Ri/Ro] 

Thus the electric field magnitude E{r), needed to calculate 
the trap cross-sections aq{r), is easily obtained. 
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IV. Trapping probability 

Trapping of a mobile, charged particle (electron or hole) 
is a stochastic process, meaning that the probability that a 
particle at r' will be trapped in the infinitesimal distance 
interval [r',r' + dr'] is pq{z)aq{r')dr' . Then the probability 
that it will not be immediately trapped is 1 — pq{z)aq{r')dr' , 
which effectively equals exp [—pq{z)aq{r')dr']. By taking the 
product of many such exponentials, the probability pq for the 
charged particle created at r, to successfully reach r is 



Pq iri,r,z) = exp 



(13) 



where use of the absolute value allows for the case (r < r^) 
that the particle moves towards the inner electrode. 

The probability for the particle, having been created at 
the interaction point r^, to be subsequently trapped in the 
infinitesimal interval [r, r + dr] is then — "^Pi^^'j^'^^ dr. To see 
this, note that 



dr 



dpq {n,r',z) , X , / \ 

-dr = -pq{ri,r,z)+pq{ri,r^,z) 



dr' 



i-Pq {n,r,z) 



(14) 



is the probability that a particle, created at r^, never arrives at 
r. Thus the derivative _^££(li£i£i ^ Xq {ri,r,z) is the PDF 
(probability distribution function) for the particle trap position 
r given r^. 

To perform a computer simulation of particle creation and 
trapping, a trap position r is randomly selected from this distri- 
bution. How is this done? The formula for converting a random 
number x taken from the uniform probability distribution 
P{x) — 1 (such X values are produced by standard random 
number generators) to the corresponding r value is derived 
as follows. The probabilities Tq {ri, r, z) dr and P{x)dx must 
be equal, so Tq {ri,r, z) dr = dx. Then integrating the former 
from ri to r, and the latter from to x gives 



X = I Tq {ri,r',z)dr' 

dpq {r,,r',z) , 

dr ^l~pq [r,, r, z) 



(15) 



which relates a randomly chosen x value to an r value. The 
function x{r) ~ 1—pq {ri,r, z) must be inverted so as to give 
r when x is chosen randomly from the interval [0, 1); that is, 
the function r{x) must be found. This is done by expressing 
the relation as 



In [1 - x] = In [pq {ri,r, z)] = - 

dr' 



Pq{z)aq{r')dr' 



= -Pq{z) 



16eo 



E{r') 



(16) 



This integral can be solved analytically by noting that 
i?(r)~'^ = _A^£_ 2 7 giving (after some careful algebra) 



e Paid 



In 



2e' 



P_„2 



(17) 



under the condition that p ^ (i.e., that the crystal is doped). 

If the X value chosen from the interval [0,1) is greater 
than 1 — Pq {ri, R, z), where R is the radius of the contact to 
which the particle is moving, then the particle has successfully 
reached that contact. That is, if the x value satisfies the relation 



\n[l~x] < ±S.p Jz)-^ 
e Paid 



In 



-A+f^R' 



-A 



2e' i 



(18) 



then the particle has successfully reached the contact at radius 
R. Otherwise the position r at which the particle is trapped is 
obtained from the equation 

= — {a+ (-^+ ^rf) (1 -a;)'^^^^;;^}. (19) 



A similar stochastic approach is taken for obtaining r^, the 
radial location at which a gamma-ray enters the crystal. For 
simplicity, the gamma-ray is assumed to shed n electron-hole 
pairs at random points (r^, z) as it traverses the length of the 
crystal. Since the z-axis of the detector points at the gamma- 
ray source, the areal distribution of in a simulation should 
be uniform over a z-slice of the crystal. Thus the areal density 
of ri points is constant; call it pi (points per area). Then dx = 
PidA giving 



2TTrdr 



Rl 



Inverting the function x [ri ) gives 

r, - {ri + x{r:i 



Rl)) 



1/2 



(20) 



(21) 



from which a value is obtained by randomly choosing a 
value X from the interval [0, 1]. 

V. Simulation algorithm 

The pieces developed above are assembled into a computer 
model of a coaxial HPGe detector. The inputs are the param- 
eter values for the detector (i?Oi Ri,L,(j) {Rn) , (j) (Ri) , Paid), 
the spectral line {E^) of interest, and the neutron fluence N/A. 
Then the model considers the gamma-rays emitted from a 
source to be normally incident on the z = (top) surface of 
the detector; subsequently each gamma-ray maintains its radial 
position ri and produces n electron-hole pairs as it traverses 
the length of the crystal. For each gamma-ray, the values 
and n are obtained stochastically according to Eq. (ISTT i and as 
described at the end of section II, respectively, and the n pair 
creation points (r^, z) are distributed randomly over the length 
L (that is, the z value for a pair is taken randomly from the 
interval [0,i] of the uniform distribution). 

The contribution of each electron-hole pair to the recorded 
energy of the gamma-ray is obtained by the following steps; 
(i) The trap density Pq{z) is calculated for each particle by 
Eq. (|9]l. (ii) This allows the terminal positions r^ and rh to 
be obtained stochastically according to Eqs. ( fTsT i and ( fT9l ). 
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1324 1326 1328 1330 1332 1334 
Gamma-ray energy (keV) 

Fig. 1. Peak coiresponding to the 1332 keV ^''Co spectral line obtained by 
the p-type detector model for a neutron fluence of 10^ cm~^. This model 
detector has the same dimensions and bias potential as the p-type detector 
studied by Pehl et al. [D. 




1260 1270 1280 1290 1300 1310 1320 1330 1340 
Gamma-ray energy (keV) 

Fig. 2. Peak corresponding to the 1332 keV ^''Co spectral Hue obtained by 
the p-type detector model for a neutron fluence of 10^'' cm~^. This model 
detector has the same dimensions and bias potential as the p-type detector 
studied by Pehl et al. (1]. 



(iii) Using these values and r/i, the induced charge Qo 
at the inner electrode is calculated by Eq. (|5]|. (iv) Then the 
contribution of the electron-hole pair to the recorded energy is 
(IQol /e) X 2.96 eV. The contributions of all n pairs constitute 
the recorded energy of the gamma-ray. 

The energy peak is constructed as a histogram of the 
gamma-ray energies. The finite width of the peak is due to 
the variation in the n value (around the average value n^) 
for gamma-rays producing the peak, and to hole and electron 
trapping suffered by some fraction of the n pairs produced by 
each gamma-ray. Thus the width of the peak can be modified 
by adjusting the values of the dimensionless parameters at 
and aq, or rather, the product at ■ aq = Aq. According to 
section III, the value Ah should lie in the range [0.1, 1], and 
Ae should be two orders of magnitude or so smaller 

This "tuning" is accomplished by reproducing the data of 
R. H. Pehl et al. |1 1. As mentioned above, two HPGe coaxial 
detectors (one n-type, the other p-type), fabricated from the 
same crystal, were irradiated simultaneously with fast neutrons 
from an unmoderated ^^^Cf source. The input to the model is 
the following: inner radius Rq = 4 mm; outer radius Ri = 21 
mm; crystal length L = 30 mm; bias potential for the p-type 
detector (i?i) = 1.6 kV; bias potential for the n-type detector 
(j) (Rq) — 2.8 kV. As dopant densities are not provided, the 
typical value pa/d — 10^" cm"'^ is used. The key data point 
from Ref. lH is the FWHM resolution of 6 keV for the 1332 
keV ^°Co line obtained by the p-type detector after a neutron 
fluence of 10^ cm~^. This experimental result is reproduced 
by the model when Ah = 0.30 (and Ag — 0.001), as shown 
in Fig. [T] 

As a check (since in fact the energy resolution is very 
sensitive to the value of Ah), the model gives the FWHM 
resolution of 64 keV after a neutron fluence of 10^° cm^^ (as 
shown in Fig. |2]i, to be compared with the experimental result 
of 70 keV. Figures |3] and |4] show the corresponding model 
results for the n-type detector (model/experimental FWHM 
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Fig. 3. Peak corresponding to the 1332 keV ^''Co spectral hne obtained by 
the n-type detector model for a neutron fluence of 10^ cm~^. This model 
detector has the same dimensions and bias potential as the n-type detector 
studied by Pehl et al. [ij. 



resolutions of 2.1/1.95 keV and 2.7/2.7 keV, respectively). 

Figures |5] and |6] show the model results for the p-type and 
n-type detectors, respectively, after a neutron fluence of 10^ 
cm"^ (model/experimental FWHM resolutions of 1.85/2.1 
keV and 1.80/1.8 keV, respectively). For this low fluence 
the model FWHM resolutions are nearly identical for the two 
detectors; instead the effect of hole trapping shows up in the 
magnitude of the shift of the peak centroid away from 1332 
keV. 

VI. Application to the INL Micro-Detective 

We exposed an ORTEC Micro-Detective (p-type) detector 
to a neutron fluence (from a ^^^Cf source) up to 10^ cm~^. 
This exercise was intended to determine the neutron fluence 



6 




U 100 



1315 1320 1325 1330 1335 
Gamma-ray energy (keV) 

Fig. 4. Peak coiresponding to the 1332 keV ^''Co spectral line obtained by 
the n-type detector model for a neutron fluence of lO^'' cm~^. This model 
detector has the same dimensions and bias potential as the n-type detector 
studied by Pehl et ai. H). 
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Fig. 6. Peak corresponding to the 1332 keV ^^Co spectral line obtained by 
the n-type detector model for a neutron fluence of 10* cm~^. This mode! 
detector has the same dimensions and bias potential as the n-type detector 
studied by Pehl et al. (TJ. 
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Fig. 5. Peak corresponding to the 1332 keV ^'^Co spectral line obtained by 
the p-type detector model for a neutron fluence of 10* cm~^. This model 
detector has the same dimensions and bias potential as the p-type detector 
studied by Pehl et al. [I]. 
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Fig. 7. Measured FWHM resolution of the 1332 keV ^"Co peak obtained 
by the Micro-Detective and GMX detectors after various neutron fluences 
(where two values are shown for the same fluence, the higher/lower value 
was measured before/after a thermal cycle). The calculated values obtained 
by the Micro-Detective model are indicated by the open circles. 



at which the energy resolution of the Micro-Detective was 
too degraded to allow its use in place of our GMX (n-type) 
detector. Figure |7] shows values of the FWHM resolution of 
the 1332 keV ^°Co peak obtained by the Micro-Detective 
and GMX detectors at various fluences (where two values 
are provided for the same fluence, the higher/lower value was 
measured before/after a thermal cycle). These results provide 
an opportunity to test the phenomenological model. 

The input to the Micro-Detective model is the following: 
inner radius i?o = 4.5 mm; outer radius i?i = 25 mm; crystal 
length L — mm; bias potential <f>{Ri) = 3.0 kV. As the 
dopant density is not provided, the typical value = 10^" 
cm~'^ is used. The calculated values of the FWHM resolution, 
indicated by the open circles in Fig. |7] compare well with 



the trend of measured values. Note that the energy resolution 
of the Micro-Detective after a neutron fluence of 10^ cm^^ 
is better than that of the p-type detector studied by Pehl et 
al. (5.4 keV versus 6 keV) despite its significantly larger size 
(diameter 50 mm versus 42 mm), due to its larger bias potential 
(3 kV versus 1.6 kV) which, by producing a stronger electric 
field across the crystal, reduces the hole trap cross-section. 

The input to the GMX (n-type) detector model is the 
following: inner radius Rq = 5.6 mm; outer radius i?i = 31.4 
mm; crystal length L = 70.3 mm; bias potential (p {R-o) — 5-0 
kV. As the dopant density is not provided, the typical value 
Pd — 10^° cm~'^ is used. The calculated FWHM resolutions 
of the 1332 keV peak at neutron fluences of lO'* cm^^ and 10^ 



cm ^ are 1.7 keV and 2.0 keV, respectively, in good agreement 
with the measured values in Fig. Q 

VII. Concluding remarks 

The main attributes of this model of neutron-damaged 
coaxial HPGe detectors are (i) the use of induced charge at 
the electrodes to determine the contribution of an electron- 
hole pair to the measured gamma-ray energy, and (ii) the use 
of stochastic methods to simulate what are, in fact, stochastic 
processes. The data of R. H. Pehl et al. [IJ provided a value 
for the dimensionless parameter Ah, related to the probability 
for immediate trapping of a mobile hole, needed to complete 
the phenomenological model. As the model is, for the most 
part, one-dimensional, it is easy to implement in a computer 
code. However, by ignoring the detector "cap" (where the 
applied electric field is not purely radial), this model is not 
well suited for application to low-energy gamma-rays which 
may be substantially stopped in that volume. 

Some observations: (i) The n-type detectors maintain good 
energy resolution to neutron fluences of at least 10^ cm^^. 
The noticeable effect of neutron damage is to shift the peak 
centroid to a lower energy, (ii) The peak shapes after a high 
neutron fluence of 10^" cm^^ are very different for p-type and 
n-type detectors. This is due not only to more hole trapping in 
a p-type detector, but also to the fact that an electron-hole pair 
with a trapped hole near the outer contact induces a smaller 
charge \Q\ at the two electrodes of a p-type detector than 
at the electrodes of a same-sized n-type detector (see section 
II for a more precise discussion of this effect). Thus those 
gamma-rays that interact closer to the outer contact (which is 
to say, most of the gamma-rays) are more likely to register 
as "counts" at the low/high end of the energy spectrum in the 
case of p-type/n-type coaxial detectors. 
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